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Abstract 

Recent theoretical and experimental studies have shown that Self- Amplified Spon- 
taneous Emission Free Electron Laser (SASE FEL) with a planar undulator holds a 
potential for generation of relatively strong coherent radiation at the third harmonic 
of the fundamental frequency. Here we present detailed study of the nonlinear har- 
monic generation in the SASE FEL obtained with time-dependent FEL simulation 
code. Using similarity techniques we present universal dependencies for temporal, 
spectral, and statistical properties of the odd harmonics of the radiation from SASE 
FEL. In particular, we derived universal formulae for radiation power of the odd 
harmonics at saturation. It was also found that coherence time at saturation falls in- 
versely proportional to harmonic number, and relative spectrum bandwidth remains 
constant with harmonic number. 
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1 Introduction 



Radiation of the electron beam in the planar undulator contains rich harmonic 
spectrum. This referes to both, incoherent and coherent radiation as well. Dur- 
ing last years a significant efforts of researchers have been devoted for studying 
the process of the higher harmonic generation in the high-gain free electron 
lasers [1]-[1 1] . Such an interest has been mainly driven by practical needs for 
prediction of the properties of X-ray free electron lasers. A fraction of a higher 
harmonic content is very important for users planning experiments at X-ray 
FEL facility. On the one hand, higher harmonics constitute rather harmful 
background for some classes of experiments. On the other hand, higher har- 
monic radiation can significantly extend operating band of the user facility. 
In both cases it is highly desirable to know properties of the hiher harmonic 
radiation. Analytical techniques have been used to predict properties of the 
higher harmonics for FEL amplifier operating in the linear mode of operation 
[8,9]. However, the most fraction of the radiation power is produced in the 
nonlinear regime, and a set of assumptions needs to be accepted in order to 
estimate saturation power of higher harmonics on the base of extrapolation of 
analytical results. As for statistical properties, they could not be extrapolated 
from linear theory at all. A lot of studies has been performed with numerical 
simulation codes. These studies developed in two directions. The first direc- 
tion is investigations of higher harmonic phenomena by means of steady-state 
codes [4,5,6,7]. Despite the results of these studies are applicable to externally 
seeded FEL amplifiers, it is relevant to appreciate that they gave the first pre- 
dictions for high radiation power in higher harmonics of SASE FEL. Another 
direction was an extraction of time structure for the beam bunching from 
time-dependent simulation code with subsequent use of analytical formulae of 
the linear theory [8]. Giving an estimate for the power, such an approach does 
not allow to describe statistical properties of the output radiation. 

In this paper we perform comprehensive study of the statistical properties 
of the odd harmonic radiation from SASE FEL. The study is performed in 
the framework of one-dimensional model with time- dependent simulation code 
FAST [12,14] upgraded for simulation of higher harmonic generation. We re- 
strict our study with odd harmonics produced in the SASE FEL. We omit from 
consideration an effect of self-consistent amplification of the higher harmonics. 
In other words, we solve only electrodynamic problem assuming that particle 
motion is governed by the fundamental harmonic. The latter approximation is 
valid when power in higher harmonics is much less than in the fundamental. 
This does not limit practical applicability of the results: it has been shown in 
earlier papers that the growth rate of higher harmonics is too small to pro- 
duce visible increase of the coherent amplification above shot noise in X-ray 
FELs [8] . Under this approximation and using similarity techniques we derive 
universal relations describing general properties of the odd harmonics in the 
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SASE FEL: power, statistical and spectral properties. The results are illus- 
trated for the 3rd and 5th harmonic having practical importance for X-ray 
FELs. 



2 Basic relations 



The one-dimensional model describes the amplification of the plane electro- 
magnetic wave by the electron beam in the undulator. When space charge 
and energy spread effects can be neglected, operation of an FEL amplifier 
is described in terms of the gain parameter T = [njoKf/^lAXwl 3 )} 1 ^ , ef- 
ficiency parameter p = A w r/(4-7r), and detuning parameter C = [27r/A w — 
lu(1 + fT 2 /2)/(2c7 2 )]/r (see,e.g. [13,14]). Here A w is undulator period, K = 
e\ w H w /2nmc 2 is undulator parameter, 7 is relativistic factor, H w is undula- 
tor field, j is the beam current density, (— e) and m are charge and mass of 
electron, I a = mc 3 /e ~ 17 kA, and u is frequency of electromagnetic wave. 
Coupling factor Kh is given by 



Q = K 2 /[2(l + K 2 )], and K is rms undulator parameter. When describing 
start-up from shot noise, one more parameters of the theory appears - number 
of particles in coherence volume, N c = I/(epcu), where I is beam current. 

We do not present here general technical details of the time dependent simu- 
lations, they have described previously in details [12,14]. The only add-on to 
these description is particle loading tool, but it is similar to that described in 
other papers (see, e.g. [15] and references therein). We note only that under 
accepted approximation (particle's dynamics is governed by the fundamental 
harmonic) we can simply calculate odd harmonics from particle distribution, 
and amplitude of the electric field scales as 



where ah is h-ih harmonic of the beam bunching. Thus, we find that coupling 
factor Kh, and time- dependent integral of the beam bunching become to be 
factorized. This allows us to extract universal ratio of the power of higher 
harmonics to the power of fundamental harmonic. 



K h = K{-l)^l 2 [J {h „ m {Q) - J (h+m (Q)} , 
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3 Statistical properties of the odd harmonics of the radiation from 
SASE FEL 

In this section we present the results of numerical studies of the operation of 
the SASE FEL in the linear and nonlinear regimes. In the framework of the ac- 
cepted model, the input parameter of the system is the number of cooperating 
electrons N c . Most of the statistical characteristics of the SASE FEL process 
are functions of iV c only in the fixed z coordinate. A typical range of the values 
of N c is 10 6 — 10 9 for the SASE FELs of wavelength range from X-ray up to 
infrared. The numerical results, presented in this section, are calculated for 
the value JV C = 3x 10 7 which is typical for a VUV FEL. It is worth mentioning 
that the dependence of the output parameters of the SASE FEL on the value 
of N c is rather weak, in fact logarithmic. Therefore, the obtained results are 
pretty general and can be used for the estimation of the parameters of actual 
devices with sufficient accuracy. 

3.1 Temporal characteristics 

Figure 1 presents a typical time structure of the 1st and the 3rd harmonic of 
the radiation from a SASE FEL at different undulator length z = Tz = 10 — 13. 
Normalized power of h-th harmonic is defined as f\ h = W h x (Ki/ K h ) 2 / (pW h ). 
Longitudinal coordinate along the pulse is s = puj (z/v z —t). The head of the 
pulse is located in the positive direction of s. A plot for the averaged power of 
the 1st harmonic is shown in Fig 2 with a solid line. It is seen that saturation 
is achieved at the undulator z = 13. Saturation length is well described in 
terms of the number of cooperating electrons N c [14,16]: 

i sat ~3 + -^=lniV c . (3) 

The normalized efficiency at saturation, i7 sat = W / sat /(pW / b) — 1.08, is almost 
independent of the value of iV c . Dashed and dotted lines show a normalized 
power ratio, fjh/vi — (Wh/Wi) x (Ki/K h ) 2 , for the 3rd and the 5th harmonic. 
One can notice that power of the higher harmonics becomes to be above the 
shot noise level only in the end of linear regime. This becomes clear if one 
takes into account that the shot noise level of the beam bunching is about 
1/y/Nc. We consider an example typical for VUV FEL with N c = 3 x 10 7 
which corresponds to the shot noise beam bunching a ~ 2 x 10~ 4 . When 
FEL amplifier operates in the linear regime, odd harmonics grow as a^, and 
we expect from this simple physical estimation that coherent contribution 
into higher harmonics can exceed the shot noise level only for the values of 
the beam bunching at the fundamental harmonic ai > 0.1, i.e. in the end 
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Fig. 1. Normalized power in the radiation pulse versus s = puo(z/v z —t) at different 
lengths of the FEL amplifier z = 10 — 13. Left and right columns correspond to the 
fundamental and 3rd harmonic, respectively 



of the linear regime. Note that shot noise level becomes to be higher when 
approaching to X-ray region. 
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Fig. 2. Normalized averaged power of a fundamental harmonic of SASE 
FEL, r)i = Pi function of a normalized undulator length 

(solid line). Dashed and dotted lines show a normalized power ratio, 
Vh/vi = (Wh/Wi) x (Ki/Kh) 2 , for the 3rd and 5th harmonic 
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Fig. 3. Ratio of coupling factors, (Kh/Ki) 2 , for the 3rd (solid line) and the 5th 
(dashed line) harmonics with respect the fundamental harmonic versus rms value 
of undulator parameter K Tms 

The plots presented in Fig. 1 allows to trace the evolution of the 3rd harmonic 
power from z — 10 (when it just started to exceed shot noise level) up to satu- 
ration point z = 13. When the beam bunching at the fundamental harmonic is 
governed by a sine function (i.e. SASE FEL operates in the linear regime), we 
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Fig. 4. Normalized power growth rate for the 1st, 3rd, and 5th harmonic (solid, 
dashed, and dotted line, respectively) 

occur well-known mechanism of the higher harmonic generation, i.e. ah oc aj - , 
and spikes of the 3rd harmonic radiation become rather pronouncing. At this 
stage, for instance, the growth rate of the higher harmonics should be propor- 
tional to the harmonic number [4,5,8]. The plot of the power growth rate for 
the 1st, 3rd, and 5th harmonic is shown in Fig. 4. One can see from this plot 
that in practical situation for SASE FEL the growth rates of higher harmon- 
ics are visibly less than those given by prediction of the linear theory. This is 
due to the fact that the value of a\ ~ 0.1 can not be considered as a linear 
stage, and the beam density modulation is not a sine-like due to nonlinear 
effects. The noise nature of the SASE FEL makes a big difference in the be- 
havior of the growth rates with respect to preidctions given in the framework 
of steady-state simulations [4,5]. Analyzing the plot for the power growth rate 
we can state that in practical situation predictions of the steady-state theory 
are valid only for the 3rd harmonic, but only on a short piece of undulator close 
to saturation, of about one gain length. Quantitative analysis of this stage of 
amplification shows that a prediction for the relation between averaged val- 
ues of the beam bunching at the third harmonic, < 1 0.3 1 2 >= 6 < |ai| 2 > 3 , 
holds approximately, and is strongly violated for higher harmonics, because of 
strong contribution of the shot noise. This feature of the SASE FEL has been 
highlighted qualitatively in early papers [8] with analysis of simulation results 
obtained with code GINGER [3]. Here we just presented more quantitative 
study. 

The plots in Fig. 2 present a general result for a ratio of the power in the 
higher harmonics with respect to the fundamental one. For the saturation we 
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Fig. 5. Normalized power ratio at saturation, (Wh/Wi) x (Ki/Kh) 2 , for the 3rd 
(solid line) and 5th (dashed line) harmonic as a function of energy spread parameter 
A^. SASE FEL opeartes at saturation 

find a universal dependency: 



^-yUt - 0.094 x-^, 



^yU-0.03x-^. 



(4) 



Universal functions for the ratio (Kh/Ki) 2 are plotted in Fig. 3. Asimptotic 
values for large value of undulator parameter are: (K^/Ki) 2 ~ 0.22, and 
(K 5 /K 1 ) 2 ~ 0.11. Thus, we can state that contribution of the 3rd harmonic 
into the total radiation power of SASE FEL at saturation could not exceed a 
level of 2%. Thus, its influence on the beam dynamics should be small. This 
result justifies a basic assumption used for derivation of a universal relation 
(4). A contribution of the 5th harmonic into the total power at saturation 
could not exceed the value of 0.3%. 

Another important topic is an impact of the electron beam quality on the non- 
linear harmonic generation process. In the framework of the one-dimensional 
theory this effect is described with the energy spread parameter A| [14] : 

A? - « A£ > ! > 



P 2 El 



where {(AE) 2 ) is the rms energy spread. Thus, result given by (4) is gener- 
alized to the case of finite energy spread with the plot presented in Fig. 5. 
We see that the energy spread in the electron beam suppresses power of the 
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higher harmonics. Within practical range of A^ this suppression can be about 
a factor of 3 for the 3rd harmonic, and about an order of magnitude for the 
5th harmonic. For practical estimations one should use an effective value of 
the energy spread describing contribution of the energy spread and emiitance 
to the longitudinal velocity spread [14]: 

((AE) 2 ) cS _ ((AE) 2 ) 2 7z V 
El El + [3 2 ' 



where 7 2 is longitudinal relativistic factor, e is beam emittance, and (3 is fo- 
cusing beta-function. The plot in Fig. 5 covers practical range of parameters 
fo X-ray FELs. The saturation length at = 0.5 is increased by a factor of 
1.5 with respect to the "cold" beam case A| = 0. 



3.2 Probability distributions 



The next step of our study is the behavior of the probability distribution of the 
instantaneous power. In Fig. 6 we show the normalized rms deviation of the 
instantaneous radiation power, a w = ((W — (W)) 2 ) 1 ^ 2 / (W) , as a function of 
the undulator length. We see that at the initial stage of SASE FEL operation 
rms deviation of the instantaneous power is equal to one for all harmonics. 
This is a consequence of start-up from the shot noise in the electron beam: sta- 
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Fig. 6. Normalized rms deviation of the fluctuations of the instanteneous radiation 
power as a function of the normalized undulator length. Solid, dashed, and dotted 
lines correspond to the fundamental, 3rd, and 5th harmonic, respectively 



8 



tistical properties of the undulator radiation and of the radiation from SASE 
FEL operating in the linear regime are governed by Gaussian statistics [14,16]. 
One of the important features of the Gaussian statistics is that the normal- 
ized rms deviation of the instantaneous radiation power is equal to the unity. 
For the fundamental harmonic statistics of the radiation becomes to be non- 
Gaussian when the amplification process enters non-linear mode [14,16]. For 
the higher harmonics non-Gaussian statistics takes place when the nonlinear 
harmonic generation starts to dominate above incoherent radiation (at z > 8 
in the present numerical example). Analytical theory of nonlinear harmonic 
generation [8] predicts the value of er w ~ 4 for the third harmonic. Analysis 
of the relevant curve in Fig. 6 shows that this prediction holds approximately 
in a short piece of the undulator length only. As we explained above, this is 
due to the fact that nonlinear harmonic generation starts to dominate above 
incoherent radiation only at the values of the beam bunching at the funda- 
mental harmonic ai ~ 0.1. However, at such a value of the beam bunching the 
modulation of the beam density already deviates from a sin-like shape due to 
nonlinear effects. 

Probability density distributions for the instantaneous power of the funda- 
mental and the 3rd harmonic are presented in Fig. 7. SASE radiation is a 
stochastic object and at a given time it is impossible to predict the amount 
of energy which flows to a detector. The initial modulation of the electron 
beam is defined by the shot noise and has a white spectrum. The high-gain 
FEL amplifier cuts and amplifies only a narrow frequency band of the initial 
spectrum Au/u <C 1. In the time domain, the temporal structure of the funda- 
mental harmonic radiation is chaotic with many random spikes, with a typical 
duration given by the inverse width of the spectrum envelope. Even without 
performing numerical simulations, we can describe some general properties of 
the fundamental harmonic of the radiation from the SASE FEL operating in 
the linear regime. Indeed, in this case we deal with Gaussian statistics. As a 
result, the probability distribution of the instantaneous radiation intensity W 
should be the negative exponential probability density distribution [14,16]: 



Here one should realize clearly that the notion of instantaneous intensity refers 
to a certain moment in time, and that the analysis must be performed over an 
ensemble of pulses. Also, the energy in the radiation pulse E should fluctuate 
in accordance with the gamma distribution [14,16]: 




(5) 




M-l 



1 



E 



) 



(E) 



cxp —M 



(E) 



(6) 



9 






A 




z = 10 




3rd harmonic 




ll INI- Ik II iNHni jniniTirTHTimf? 



z = 1 1 

3rd harmonic 




Fig. 7. Probability distribution of instanteneous radiation power at different lengths 
of the FEL amplifier z = 10 — 13. Left and right columns correspond to the funda- 
mental and 3rd harmonic, respectively. Solid line shows probability density function 
(7) 

where T(M) is the gamma function of argument M, andl/M= ((E — (E)) 2 } / (E) 2 
is the normalized dispersion of the energy distribution. These properties are 
well known in statistical optics as properties of completely chaotic polarized 
radiation [18]. 



The statistics of the high-harmonic radiation from the SASE FEL changes 
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significantly with respect to the fundamental harmonic (e.g., with respect to 
Gaussian statistics). It is interesting in our case to be able to determine the 
probability density function of instantaneous intensity of SASE radiation after 
it has been subjected to nonlinear transformation. We know the probability 
density function p(W) = (W)^ 1 exp(— W/{W)) of the fundamental intensity 
W, and W is subjected to a transformation z = (W) n . The problem is then 
to find the probability density function p(z). It can be readily shown that this 
probabilty distribution is [17]: 

p( z ) = -J—z^-^expi-z^/iW)) . (7) 



Using this distribution we get the expression for the mean value: (z) = n\ (W) n . 
Thus, the nth-harmonic radiation for the SASE FEL has an intensity level 
roughly n\ times larger than the corresponding steady-state case, but with 
more shot-to-shot fluctuations compared to the fundamental [8]. Nontrivial 
behavior of the intensity of the high harmonic reflects the complicated nonlin- 
ear transformation of the fundamental harmonic statistics. One can see that 
Gaussian statistics is no longer valid. Upper plots in Fig. 7 give an illustration 
to these consideration. Despite in practical example we do not have pure linear 
amplification regime, the probability density functions for the instanteneous 
power follow rather well prediction (7). 

Analysis of the probability distributions in Fig. 7 shows that in the nonlinear 
regime, near the saturation point, the distributions change significantly with 
respect to the linear regime for both, the fundamental and the 3rd harmonic. 
An important message is that at the saturation point the 3rd harmonic radia- 
tion exhibits much more noisy behaviour (nearly negative exponential) while 
stabilization of the fluctuations of the fundamental harmonics takes place. 



3. 3 Correlation functions 



The first and the second order time correlation functions are defined follows: 



9i(t-t') 
92(t-t') 



(E(t)E*(t) 



(mw)(\E(t'W) 
(\E(t)\ 2 \E(m 

(\E(t)\ 2 )(\E(t'W) 



1/2 



(8) 



In Fig. 8 we show the evolution of the time correlation functions of first and 
second order. At each normalized position along the undulator, z, they are 
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Fig. 8. First (left column) and second (right column) order correlation function at 
different lengths of the FEL amplifier z = 10 — 13. Solid, dashed, and dotted lines 
correspond to the fundamental, 3rd and 5th harmonic, respectively 

plotted versus the normalized variable f = pu (t — t'). Upper plot in Fig. 8 
corresponds to the linear stage of SASE FEL operation. In the case of the 
fundamental harmonice we deal with a Gaussian random process and the 
relation between the correlation functions holds for g2{t — t') = 1 + \gi(t — 1')\ 2 . 
This feature does not hold place for figher harmonics. The nontrivial behavior 
of the correlation functions reflects the complicated nonlinear evolution of the 
SASE FEL process. The second-order correlation function of zero argument, 
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Fig. 9. Normalized coherence time of a SASE FEL as a function of normalized 
undulator length. Solid, dashed, and dotted lines correspond to the fundamental, 
3rd, and 5th harmonic, respectively 

(72(0), takes values smaller or larger than two, but always larger than unity. 
Note that there is a simple relation between #2(0) and the normalized rms 
power deviation: g 2 (0) = 1 + a\ (see Fig. 6). It is a well-known result of 
statistical optics that the cases of #2(0) = 1 and #2(0) = 2 correspond to 
stabilized single-mode laser radiation and to completely chaotic radiation from 
a thermal source, respectively. The values of #2(0) between 1 and 2 belong 
to some intermediate situation. In classical optics, a radiation source with 
#2(0) < 1 cannot exist but the case of #2(0) > 2 is possible. As one can 
see from Fig. 8, the latter phenomenon (known as superbunching) occurs for 
higher harmonics of SASE FEL, or for fundamental one when the SASE FEL 
operating in the nonlinear regime. 

In Fig. 9 we present the dependence on the undulator length of the normalized 
coherence time f c = pcuoT c , where r c is 



For the fundamental harmonic the coherence time achieves its maximal value 
near the saturation point and then decreases drastically. The maximal value 
of f c depends on the saturation length and, therefore, on the value of the 
parameter N c . With logarithmic accuracy we haave the following expression 
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for the coherence time of the fundamental harmonic: 

\'c)max — Y -j^g 

One can find from Fig. 8 that coherence time at saturation for higher harmon- 
ics approximately falls inversely proprtional to the harmonic number h. 

3.3.1 Spectral Characteristics 

When comparing radiation spectra, it is convenient to use the normalized 
spectral density, h(C), defined as 

oo 

J dCh(C) = l. 

— oo 



The frequency deviation, Auj, from the nominal value of uih can be recalculated 
as Auj = —Ipuj^C. Since we consider the model of a long rectangular bunch, 
the function h(C) can be treated as the normalized spectral density of both 





Fig. 10. Normalized spectrum at different length of the undulator: z = 10—13. Solid, 
dashed, and dotted lines correspond to the fundamental, 3rd and 5th harmonic, 
respectively 
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the radiation energy and the power. 



Normalized envelope of the radiation spectrum and the first order time corre- 
lation function are connected by the relation [18]: 



The temporal structures of the radiation pulses (see Fig. 1) are used for calcu- 
lating the first order time correlation function (see Fig 8). Then the radiation 
spectra are reconstructed by Fourier transformation of the first order time 
correlation function. Figure 10 shows evolution of the radiation spectra of the 
SASE FEL radiation from the end of the linear regime to saturation. Note that 
spectrum width of the higher harmonics from SASE FEL differs significantly 
from that of incoherent radiation. For the case of incoherent radiation relative 
spectrum width, Auj/ujh scales inversely proportional to the harmonic num- 
ber h [19]. One can see that situation changes dramatically for the case when 
nonlenear harmonic genaration process starts to be dominant. At saturation 
we find that relative spectrum bandwidth becomes to be nearly the same for 
all odd harmonics. 



4 Summary 

In this paper we performed detailed study of the properties of the odd har- 
monic of the radiation from SASE FEL. Universal formulae for contribution 
of the higher odd harmonics to the FEL power for SASE FEL operating at 
saturation are obtained. In the case of cold electron beam these contributions 
are functions of the undulator parameter K only. General statistical properties 
of the odd harmonics of the SASE FEL operating in saturation are as follows. 
Power of higher harmonics is subjected to larger fluctuations than that of 
the fundamental one. Probability distributions of the instantaneous power of 
higher harmonics is close to the negative exponential distribution. The coher- 
ence time at saturation falls inversely proportional to harmonic number, and 
relative spectrum bandwidth remains constant with harmonic number. 
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